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Muon-spin rotation study of the in-plane magnetic penetration depth of FeSeo.85: 

evidence for nodeless superconductivity 
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The in-plane magnetic penetration depth Aab of the iron selenide superconductor with the nominal 
composition FeSeo.85 was studied by means of muon-spin rotation. The measurements of \~^{T) are 
inconsistent with the presence of nodes in the gap as well as with a simple isotropic s— wave type of 
the order parameter, but can be equally well described within a two-gap (s -I- s) and an anisotropic 
s— wave gap picture. This implies that the superconducting energy gap in FeSeo.gs contains no 
nodes. 

PACS numbers: 76.75.-|-i, 74.70.-b 



The recent discovery of the Fe-based high-temperature 
superconductors has attracted considerable attention to 
the pnictides. The superconductivity was first found 
in LaOi_a:Fa:FeAs [l| and, later on, in other single- 
layer arseno-pnictides by replacing La with various rare- 
earth elements (Sm, Nd, Pr, Ce, Ho, Y, Dy, and Tb) 
as well as in oxygen-free compounds such as doped 
double-layer MFe2As2 (M=Ba, Sr, and Ca) 3 and single- 
layer LiFeAs [1]. The common structural feature of all 
these materials is the Fe-As layers consisting of an Fe 
square planar sheet tetrahedrally coordinated by As. Re- 
cently, superconductivity with the transition tempera- 
ture Tc ~ 8 K was discovered in a— FeSe with PbO- 
structure Q . This compound also has a Fe square lattice 
with Fe atoms tetrahedrally coordinated by Se ones, sim- 
ilar to the structure of FeAs planes in the single- and the 
double-layer arseno-pnictides. In this respect FeSc, con- 
sisting of the "superconducting" Fc-Se layers only, can 
be treated as a prototype of the known families of Fe- 
As based high-temperature superconductors and, conse- 
quently, becomes a good modeling system to study the 
mechanism leading to the occurrence of superconductiv- 
ity in these new class of materials. 

Here we report a study of the in-plane magnetic 
field penetration depth Xab in iron selenide supercon- 
ductor with the nominal composition FeSeo.ss by means 
of muon-spin rotation (/iSR). X~b{T) was reconstructed 
from the temperature dependences of the /iSR lincwidth 
measured in a magnetic field of 0.01 T. The observed 
A~jf(T) was found to be equally well described within 
the framework of anisotropic s— wave as well two-gap 
s -\- s— wave models. In of anisotropic s— wave 

model the maximum value of the gap at T = was 
found to be Aq = 1.35 meV leading to 2Ao/fcs7c = 3.79, 
close to the weak-coupling BCS value 3.52. The two- 
gap s + s— wave model yields Ao,i = 1.63 meV and 



Ao,2 = 0.38 meV. The corresponding gap to Tc ratios 
are 2Ao,i/fci3rc = 4.59 and 2Ao,2/fcBTc 1.07 close 
to those reported for various single- and double-layer 
arseno-pnictide superconductors [fll3, S] • 

Details on the sample preparation for FeSeo.ss can 
be found elsewhere [51]. X-ray diffraction analysis re- 
veals that the a-FeSe phase is dominant and that the 
amount of the impurity fraction docs not exceed ~7-10%. 
The AC magnetization (Mac) measurements (poHac = 
0.1 mT, ly = 1000 Hz) were performed on a Quantum De- 
sign PPMS magnetometer at temperatures ranging from 
1.75 K to 300 K. The superconducting transition tem- 
perature Tc = 8.26(2) K was obtained as an intersection 
of the linearly extrapolated Mac (T) with Mac — const 
line (see Fig. [T]) . 




FIG. 1: (Color online) Temperature dependence of the AC 
magnetization Mac {i-loHac = 0.1 mT, 1/ = 1000 Hz) of 
FeSeo.85- 

Zero field (ZF), longitudinal field (LF) and trans- 
verse field (TF) fiSK experiments were performed at the 
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7rM3 beam line at the Paul Scherrer Institute (Villigen, 
Switzerland). Here LF and TF denote the cases when 
the magnetic field is applied parallel and perpendicular 
to the initial muon-spin polarization, respectively. The 
experiments down to T ~ 1.5 K were performed on the 
low-background General Purpose Surface-Muon instru- 
ment (GPS) and those down to T ~ 0.02 K on the Low- 
Temperature-Facility (LTF) instrument. 



This proves that the magnetism in FeSeo.ss sample stud- 
ied here is static in origin and is caused by diluted and 
randomly distributed magnetic moments. Bearing this 
in mind and by taking into account the presence of the 
relatively high paramagnetic contribution at T > Tc (see 
Fig. [T|) we may conclude that the magnetism observed in 
both, ZF /iSR and magnetization experiments, has sim- 
ilar sources and, most probably, caused by the traces of 
Fe impurities 
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FIG. 2: (Color online) Temperature dependence of the ZF 
muon depolarization rate A of FeSeo.85. The inset shows ZF 
(T = 1.5 K and 10 K) and LF (T = 10 K, fioH = 0.01 T) 
/.tSR time-spectra of FeSeo.ss. 

First we are going to present the results of the ZF ex- 
periments. In the whole temperature region the ZF data 
were found to be well described by the single-exponential 
decay function: 

A^^{t) = Ancxp{-At). (1) 

Here Aq is the initial asymmetry at t = and A is the cx- 
ponentional depolarization rate. The results of the anal- 
ysis of the ZF data and the representative ZF and LF 
muon-time spectra are shown in Fig. O The open and 
the closed symbols are from the measurements taken on 
the LTF and the GPS instruments, respectively. 

The exponential character of the muon polarization 
decay might be explained either by existence of fast 
electronic fluctuations measurable within the /xSR time- 
window or by a static magnetic field distribution 
caused by diluted and randomly oriented magnetic mo- 
ments [10|. To distinguish between these two cases LF 
/iSR experiments were performed. As is shown in Ref. Ill 
in a case when the applied longitudinal field is much 
stronger than the internal field {B > lOBint) the muon 
spins become "decoupled" from the static internal field. 
On the other hand, field fluctuations perpendicular to the 
applied external field can cause irreversible spin-fliptran- 
sitions of the muon spin, leading to depolarization [9| . As 
is shown in the inset of Fig.[2]an external field of 0.01 T is 
already enough to completely decouple the muon spins. 
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FIG. 3: (Color online) Temperature dependence of the Gaus- 
sian depolarization rate a at jj,oH = 0.01 T of FeSeo.ss. The 
inset shows the TF muon-time spectra above (T — 10 K) and 
below {T = 1.5 K) the superconduting transition temperature 
Tc = 8.26 K. 

The in-plane magnetic penetration depth Xat was stud- 
ied in the TF /iSR experiments. In a powder sample 
the magnetic penetration depth A can be extracted from 
the Gaussian muon-spin depolarization rate asc{T) oc 
1/A^(T), which probes the second moment of the mag- 
netic field distribution in the mixed state [13] . For highly 
anisotropic layered superconductors (like the pnictide su- 
perconductors) A is mainly determinated by the in-plane 
penetration depth Aab (Jsc{T) oc l/\l^,{T). By tak- 
ing into account the weak magnetism observed in our 
ZF experiments (see Fig.[2]and discussion above) the TF 
/iSR data were analyzed by using the following functional 
form: 

A^^{t) = Aq exp(-At) exp(-(T2t2/2) cos{"fB,„tt + 0). 

(2) 

Here 7/27r = 135.5 MHz/T is the muon gyromagnetic 
ratio, (f) is the initial phase of the muon-spin ensemble, 
and (T = (tTgc -|- (Tnm,)^'^ is the Gaussian relaxation rate. 
anrn IS the nuclear magnetic dipolar contribution which is 
generally temperature independent [l3|. The whole set 
of 0.01 T TF ^SR data was fitted simultaneously with 
Aq, A, and (/> as a common parameters and a and Bint as 
individual parameters for each temperature point. The 
exponential relaxation rate was assumed to be temper- 
ature independent in accordance with the results of our 
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FIG. 4: (Color online) Temperature dependence of of FeSeo.85 obtained from measured CTsc by means of Eq. ([Sjl. The 
fitting curves (solid black lines) were obtained within the following models of the gap symmetries: s— wave (a), s + s— wave (b), 
anisotropic s— wave (c), and d— wave (d). The corresponding angular dependences of the gaps are shown as the insets. 



ZF /xSR experiments (see Fig. [2]). anm was fixed to the 
value obtained above Tc where a = anm (see Fig. [3]). The 
results of the analysis and the representative TF muon- 
time spectra are shown in Fig. [31 The open and the closed 
symbols are again from the measurements taken on the 
LTF and the GPS instruments, respectively. 

The superconducting part of the Gaussian depolariza- 
tion rate asc can be converted into Xab via 0, [ij, : 



(3) 



where $o = 2.068 • 10~^^ Wb is the magnetic flux quan- 
tum. Fig. m shows A~{f(T) obtained from the mea- 
sured asc{T) by means of Eq. Regarding the pair- 
ing symmetry, available experimental results on various 
single- and double-layer arseno-pnictides are divided be- 
tween those favoring an isotropic [3, 17 1 as well as an 
anisotropic [ist nodeless gap and those supporting line 
nodes 19|. The two-gap behavior was also reported in 
Refs. Hn^, [l^, [2^, [U. Bearing this in mind the data 
in Fig. [4] were analyzed by using single-gap and two-gap 
models, assuming that the superconducting energy gaps 
have the following symmetries: s— wave (a), s + s— wave 
(b), anisotropic s— wave (c), and d— wave (d). 

Temperature dependence of the magnetic penetration 
depth A was calculated within the local (London) approx- 
imation (A ^ ^) by using the following functional form 
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(4) 

Here A ^(0) is the zero-temperature value of the mag- 
netic penetration depth, / = [1 -|- exp(i?/fcBT)]~^ is the 
Fermi function, Lp is the angle along the Fermi surface, 
and A(T, 1^9) = 1^q5{T /Tc)g{ip) (Aq is the maximum gap 
value at T = 0). The temperature dependence of the gap 
is approximated by 5{T/Tc) = tanh{1.82[1.018(rc/T - 
[231 . The function g{(p) describes the angular 
dependence of the gap and is given by g'^i^p) = 1 for 
the s— wave gap, g'^i'p) ~ \ cos{2(p)\ for the rf— wave gap, 
and g'^'^^if) = (1 -|- acos4(^)/(l + a) for the anisotropic 
s— wave gap psj . 

The two-gap calculations were performed within the 
framework of the so called a— model assuming that the 
total superfluid density is a sum of two components 
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\-\T) A-2(T,Ao,i) A-2(T,Ao,2) 

= ■ ■ r + (l-^)- - . r. (5) 



A-2(0) 



A-2(0,Ao,i) 
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Here Ao,i and Ao,2 are the zero-temperature values of the 
large and the small gap, respectively, and w (0 < u; < 
1) is the weighting factor which represents the relative 
contribution of the larger gap to A^^. 

The results of the analysis are presented in Fig. d] by 
solid black lines. The angular dependences of the gaps 
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[Aq -gifp)] are shown in the corresponding insets. It is ob- 
vious that the simple s— and d— wave approaches cannot 
describe the observed A~jf(T) [see Figs. H] (a) and (d)]. 
In both cases the low temperature points stay systemat- 
ically higher than the theoretically derived curves. The 
constant, within our experimental uncertainty, ^~^{T) 
at r < 1 K is also inconsistent with the presence of 
any type of nodes in the energy gap of FeSeo.ss- In 
contrast, both, anisotropic s— and two-gap s + s— wave 
models [Figs. |4](b), (c)] describe the experimental data 
reasonably well. In the following we are going to discuss 
separately the results obtained within the framework of 
these two models. For the anisotropic s— wave case we 
get Ao = 1.35 meV, a = 0.796, and Aab(O) = 406 nm. 
The corresponding gap to Tc ratio is 2A/kBTc = 3.79 
which is rather close to the weak-coupling BCS value 
3.52. The obtained variation with angle A™''^/A™™ = 
1.35/0.153 ~ 8.8 is substantially bigger than 1.2 re- 
ported in Ref. [H for NdFeAsOo.gFo.i. The fit within 
the two-gap s + s— wave model yields Ao,i = 1.63 meV, 
Ao,2 = 0.38 meV, w = 0.654, and Aafc(O) = 403 nm. It is 
interesting to note that the "large" and the "small" gap 
to Tc ratios 2Ao,i/fcsTc = 4.59 and 2Ao^2/kBTc = 1.07 
are very close to those reported for various single- and 
double-layer arseno-pnictide superconductors based on 
the results of the point contact Andrcev reflection spec- 
troscopy experiments of Szabo et al. j^] and Gonelli et 
al. 7 1 , and the first critical field measurements of Ren et 
al. 8| . The multiple gaps may originate from the multi- 
ple bands at the Fermi level of FeSeo.ss- First-principle 
calculation indicates that the Fermi surface (FS) of FeSe 
is quasi-two dimensional and consists of hole-type sheets 
around the F point and electron-type sheets around the 
M point of the Brillouin zone [2^. It is conceivable that 
the two gaps open up on the different sheets of the FS. 
In this context, the present compound may resemble the 
situation of MgB2 where a large gap opens on the FS de- 
rived from the orbitals in the boron plane, while a small 
gap opens on the FS derived from orbitals perpendicular 
to the boron plane f27|. 

To conclude, muon-spin rotation measurements were 
performed on the superconductor FeScg.ss {Tc — 8.3 K). 
A~jf (T) was reconstructed from the temperature depen- 
dence of the Gaussian muon depolarization rate mea- 
sured at fioH = 0.01 T. The absolute value of the in- 
plane magnetic penetration depth Xat at T = was es- 
timated to be Aa6(0) 2± 405 nm. The temperature de- 
pendence of X^^ was found to be inconsistent with an 
isotropic s— wave as well as with a wave symmetry 
of the superconducting energy gap. A good agreement 
between the experimental data and the theory was ob- 
tained within the framework of a two-gap s + s and an 
anisotropic s— wave gap models thus suggesting that the 
superconducting energy gap in FeSeg.ss superconductor 
is fully developed and contains no nodes. 
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